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ABSTRACT 



A nonlinear clutch type servomechanism is proposed to position 
the control surfaces of an underwater vehicle in accordance with specifi 
performance criteria. 

Two system configurations are discussed and compared: 

(1) A two-state system employing pulse width angular velocity 
modulation on counter rotating clutches. 

(2) A three-state system employing counter rotating clutches and 
a brake with error sampling. 

Steady state load vibration will be present in each case. The 
degrading effects of these vibrations are developed and discussed. 

Clutch power dissipation phenomena are discussed. 

Actuator specifications are examined in light of overall system 
dynamics . 

The present state of clutch technology is examined to fulfill system 
requirements . 



Thesis Supervisor: George C. Newton, Jr. 



Title: Professor of Electrical Engineering 



TABLE OF CONTENTS 

page 

CHAPTER I INTRODUCTION AND SUMMARY 1 

1.0 Problem Statement 1 

1. 1 System Configuration 1 

1.2 Actuator Configuration 3 

1.3 Control Surface Wobble 5 

1.4 Clutch Power Dissipation 7 

1 . 5 Overall System Dynamics 9 

1.6 Clutch Specifications 11 

CHAPTER II SOUND GENERATION 13 

CHAPTER HI HYDRODYNAMICS AND CONTROL 

SURFACE OSCILLATION 16 

CHAPTER IV CLUTCH POWER DISSIPATION 24 

BIBLIOGRAPHY 28 



LIST OF FIGURES 



1.1(a) Yaw-Rate Loop page 2 

1. 1(b) Control Surface Actuator 2 

1. 2(a) Two-State Clutch Servo with 

Pulse Width Modulation 8 

1.2(b) Three-State Clutch Servo Configuration 8 

1.5 Nichols Chart of Yaw-Rate 

Loop Frequency Response 10 

1.6 Comparative Clutch Characteristics 11 

3.1 Theodorsen's Function 19 

3.2 Source SPL vs. Switching Frequency 23 

4. 1(a) Analog Simulation of Clutch Dissipation 

Phenomena 25 

4. 1(b) Response of Second Order Clutch 

Output System 



26 



CHAPTER I 



INTRODUCTION AND SUMMARY 
1 . 0 PROBLEM STATEMENT 

Switching amplifiers, in general, offer higher efficiency and 
higher power per unit weight and volume than do continuous amplifiers. 
This research investigated the application of clutch type actuators 
operated in a switching mode to the feedback position control of the 
control surfaces of an underwater vehicle. The vehicle was assumed 
to be operating at maximum speed. 

1. 1 SYSTEM CONFIGURATION 

A self-propelled underwater vehicle when properly designed with 
neutral stability in yaw-rate will exhibit fastest response in that mode 
of motion. A control surface position servo satisfying the requirements 
of yaw-rate tracking will be adequate for roll and depth control. 

A linear small signal model for the yaw-rate loop is indicated 
in Fig. 1. 1(a). 

The individual transfer functions are 
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Fig. 1.1 (a) Yaw-Rate Loop 
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Fig. 1.1 (b) Control Surface Actuator 
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£, = damping ratio 

4j = command yaw rate 

6. = command rudder deflection 

i 

4j = vehicle yaw rate 

6 = actual rudder deflection 

o 

A linear model for the control surface position control system is 
indicated in Fig. 1. 1(b). The control surface system frequency response 
has been specified as flat to 7 cps with -3dB point at 27 cps.^^ The 
torque amplifier must deliver up to 1000 in -lb torque at a 

maximum speed of 100 RPM. 

1.2 ACTUATOR CONFIGURATION 

The design of a linear control system for control surface positioning 
is straightforward and can be realized with an electric or hydraulic 
actuator. These devices have inherent disadvantages, notably energy 
conversion in the former and shelf life in the latter. 

Clutch type actuators offer superior performance in many respects 
although analysis is difficult due to nonlinear describing equations. 
Counter rotating clutch actuators can be operated in a push-pull fashion 
or in a switching mode. The former is undesirable from efficiency 
and heat sinking considerations. Switching operation is considered in 
this research. 

Two basic configurations are considered: A two-state pulse width 

modulated counter-rotating clutch arrangement and a three-state 
system utilizing two counter-rotating clutches and a brake with state 
determined by sampled errors. 
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The two-state system is functionally diagrammed in Fig . 1.2(a). 

The pulse width modulator output is designed to operate at a 

constant recurrence rate with fixed output magnitude positive or negative 

(3) 

with half cycle duration determined by input level. The pulse width 
modulator output drives the actuation mechanism of the clutch. 

High inertia clutch input members are driven at constant speed 
from the vehicle prime mover. Assume that input and output members 
become instantly locked on signal from the pulse width modulator. 

The clutch output member velocity will be a square wave. The spring 
Kj. transforms this velocity square wave into a triangular torque wave 
whose short-time mean level is transmitted to the load. The load and 
filter inertia attenuate torque components at switching frequency 
and higher harmonics. 

Compensation and feedback networks insure reset action and 
stabilize the system. Reduction gearing is necessary to reduce clutch 
size. 

In the steady state the clutch output velocity will be a symmetrical 
square wave. Some wobble will be inherent in the control surface load 
at switching frequency and its harmonics. 

The three state system operation is s elf explanatory in Fig . 1.2(b). 
System operation is quite similar to that of the two-state configuration. 
The notable exceptions are the more complicated decision box and lack 
of steady state output oscillations at the expense of steady state errors, e. 
Limit cycle oscillations can occur at sampling frequency if the forward 
loop gain is high enough. 
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Detailed design of the clutch actuator system is the subject of 
a separate thesis and will not be dealt with here. 

The advantages of clutch type actuators derive from the high 
torque squared-to-inertia ratios and low control power requirements. 
Shelf life is no problem with clutch servos. 

The disadvantages of both configurations include heat dissipation 
at contact surfaces and output wobble with the possibility of exciting 
structural resonance. 

Selection of the optimum switching of sampling frequency will 
involve tradeoffs between output wobble amplitude, which will decrease 
with switching frequency, and power dissipation in the clutches, which 
will increase with switching frequency. The fundamental switching 
frequency must lie above load resonance and away from structural 
resonances . 

1 . 3 CONTROL SURFACE 

The adverse effects of control surface wobble were investigated 
and acceptability criteria were assigned. 

Control surface wobble manifests itself as a perturbation on 
the control surface angle of attack. The hydrodynamics of the situation 
were examined. Since mean incidence remains unchanged the dynamic 
lift of the control surface is not affected. Viscous drag effects were 
second order in magnitude, in fact, with the proper phasing and correct 
relationship between frequency, vehicle velocity and control surface 
chord length, propulsion will result. Cavitation inception was considered 
unlikely since the disturbance velocities were insignificant compared 
with quasi-steady levels. 
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Sound generation by the oscillating control surface proved 
to be a significant result. The sonic intensity level can be 
determined from the equation 



I(r) = 



,dL.2 Z 
(■dt ) COS 0 
' 2 3 T 

16 tt p c r 
r o o 



watts/ nr/ 



where 



L = lift force exerted by hydrofoil 

0 ~ angle between direction of observation 
and lift force 



= density of undisturbed fluid 

c = speed of sound in undisturbed fluid 
o 

r = distance from source to observer 



This relationship is developed in Chapter Two. Thus, knowing 

lift as a function of angle of attack and angle of attack as a function 

of time, the sonic intensity can be obtained. The sound pressure level (SPL) 

can be obtained from intensity. A convenient measure of sound is the 

source level (far field radiation extrapolated to r = 1M) expressed in 

-5 2 

decibel ratio to 0.002 microbar (2 x 10 Newton/Meter ). 

For a given maximum control surface rate specification (100 RPM 
in this case) the source SPL can be determined as a function of switching 
frequency in the two-state configuration. A limit cycle amplitude 
determination must be made for the three-state configuration. 

The two-state source SPL was observed to decrease 20 dB per 



decade of switching frequency. 
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1.4 CLUTCH POWER DISSIPATION 

Referring to Fig. 1.2(a) if were zero, the amount of energy 

necessary to accelerate the clutch output member from full clockwise 
angular velocity is 

= J q £1* joules 

= clutch output inertia 
Q = angular velocity of input members 

Since this operation is performed twice in a switching period, 

P, = 2J Q Z f watts 

d c sw 

f = switching freq (cps) 

The dissipation rating of the clutch must lie above this level. 

The situation can be alleviated somewhat by making as stiff as 

possible without exceeding clutch pullout torque in spring windup. An 

analog simulation bears this out (see Fig. 4. 1(b) ). This result is 

consistent with two-state system design objectives in that clutch 

pullout torque must be higher than the greatest demanded. 

In the analog simulation the filter inertia (Fig. 1.2(a)) was 

considered infinite. In the actual system J r > 5J for efficient 

J f c 

filtering so this assumption was not a bad one. The simulation is 
detailed in Chapter Four. 

Since power dissipation is a linear function of switching 
frequency and a quadratic function of input member speed the gear ratio 
should be selected to minimize input member speed. Unfortunately 
load torque requirements determine the gear ratio since the clutch 
dimensions must be kept small. Pullout torque is determined by clutch size. 
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Fig* 1.2(b) Three State Clutch Servo Configuration 
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If the control surface can be pivoted nearer to its center of 
pressure the output torque required (1000 in-lbs) will be reduced 
and the necessary reduction ratio can be reduced permitting lower 
clutch input member speed while maintaining maximum control surface 
position rate at the specified 100 RPM. 

1 . 5 OVERALL SYSTEM DYNAMICS 

Since the control surface servo is embedded within the overall 
vehicle dynamic control loop, a study was carried out to determine 
the effect of actuator bandwidth reduction on the dynamics of the 
total system. 

The vehicle exhibits fastest behavior in the yaw rate loop (Fig. 1. 1(a)). 
The control surface servo system was modelled as a second order 
system with damping ratio equal to . 

0 CO 

0 n 

" 6 . = 7 ^ 

1 S t CO S "r (O 

n n 

It can be shown that damping ratio =0.5 minimizes integral square 
error for a linear position servo for a step input. This behavior can 
be achieved by proper selection of forward loop gain and output rate 
feedback. 

Bode plots of the open yaw-rate loop were constructed for the 
specified servo and three slower servos. These data were transferred 
to a Nichols chart for comparison (Fig. 1.5). 
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Fig. 1.5 Nichols Chart of Yaw-Rate Loop Frequency Response 



It was concluded that the response speed of the control surface 
actuator could be halved without degrading the closed yaw-rate 
loop performance. 

1 . 6 CLUTCH SPECIFICATIONS 

The state of clutch technology was scanned and unfortunately 
present production line clutches do not fulfill system requirements. 
Switching rates in the range 7 5 to 150 cps are necessary for the two 
state system. This gives rise to an input-output rise time of one or 
two milliseconds. This demand has been met by very small instrument 
magnetic particle clutches. Figure 1.6 compares characteristics of three 
production devices with proper torque and speed characteristics. 



Maximum 
Clutch Torque 

Type (in-lb) 


Maximum 

Speed 

(RPM) 


Maximum 

Dimension 

(in) 


Rise 
Time 
(±n sec) 


Output 

Inertia 

(in-lb sec^xlO^) 


Dissipatio 
Rating 
i (watts) 


Solenoid- 

Operated 

Mech. 

F riction 


7500 


3.7 


30 


6. 93 


136 


Magnetic 
Particle 35 


3400 


3. 1 


30 


2.37 


75 


Solenoid- 
Actuated 
Vacuum- 60 

Operated 
Mechanical 


1800 


5 


8 


9. 3 


127 



Fig. 1.6 Comparative Clutch Characteristics 
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In addition, the magnetic particle clutch exhibits 0.6 in-lb drag 
torque at zero excitation and the vacuum -operated clutch requires a 
small vacuum pump. 

It may be concluded that the two-state actuator configuration is 
not technically feasible at present. Since the three-state system has 
not been fully investigated, a working model may be constructed using 
the present clutch technology. High speed magnetic tape drives have 
indicated a need for fast mechanical clutches and current research may 
evolve a device capable of meeting the requirements of the two-state 



actuator system. 



CHAPTER II 



SOUND GENERATION 



The intensity of sound resulting from a lifting body in fluid flow 

(4) 

has been investigated. Rayleigh has shown that any frequencies in 
a fluid flow are equal to those of the sound generated. 

Lighthill has considered the equations of motion of fluids and 
compared them with the equations of sound propagation in a medium 
at rest. The exact equations for the former may be written in the form 



a 2 P 



where 



at 

T u 

P 

P ij 

C 



v. 

i 



- C 2 V 2 

o v p 



9x. ax. 

1 J 



T .. 

IJ 



( 2 . 1 ) 



= p v.v. + p. . - C 6. . 
i J ij o lj 

= density 

= compressive stress tensor 

= sonic velocity in fluid at rest 

= velocity components direction x. (i = 1, 2, 3) 



Physically, this means that sound is generated by a fluid flow 
exactly as in a uniform medium at rest which is acted upon by externally 
applied fluctuating stresses. If we regard the flow as known it is possible 
to evaluate the sound generated. 

Equation (2. 1) is the non-homogeneous wave equation. The most 
general solution (Stratton 1941) is 



p(x ) - p o = 




a 2 T, 



8y. ay 



ij d y 



i 'j i x_ y 



± + K A £ I? dS(y) 



r 0n 



an 



c r 0n 8t 
o 
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In this equation 



JLLii JLe_ p !_£. 

8y.ay. ’ 3n ’ ^ ’ at 

i J 



are taken at retarded times 



t - — , where r = |x - y | and n is the outward normal from the 
o 

fluid. The first integral is taken over the total volume V external 
to the solid boundaries and the second integral is taken over the 
surface S of the solid boundaries. 

Curie has shown that the solution may be simplified under the 
assumption that: 

(1) | x*] » \ ("far field" approximation) 

(Z) the solid boundaries are rigid, i.e., non - dilatational 
(3) the flow is incompressible (Mach No. << 1) 

Under these assumptions the solution reduces to 



p - p G = p = 



i 



A JL f (t) 

3 2 at i l 1 



4 ttc x 
o 



where 



V*> = / • 



t) d S(y ), the total resultant force 



exerted on the fluid by the solid boundaries. 

In acoustic media, the sound pressure field is proportional to 
the density fluctuation field by the constant c^ thus 



P = 



i x. 

i i a 



4 TTC Z at 1 

O X 



F.(t) 



This is the radiation field of a dipole source with axis parallel 
to F(t) the net resultant force exerted on the fluid by the boundaries, 
Rewriting in more convenient variables : 



A 

P 
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1 

4 TT C 
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d_ 

3t 



F.(t) 



(^. 2 ) 



o 

where 0 is the angular distance between the radius vector r and the 
axis of F.(t). 

i 

Since intensity is defined as the acoustic power transmitted through 
a unit area 



P c 
r o o 



p = mean square acoustic pressure (Newtons/M ) 



P c 
r o o 



characteristic impedance of acoustic medium 
(MKS RAYLS) 



Generally it is more physically meaningful to consider a decibel 
ratio referred to a datum level for intensity and pressure measurements. 

A customary quantitative measurement in underwater acoustics is 
the source level sound pressure level. This is defined as 

~ 2 

SPL = 10 log 1 0 -Eg- dB ' 

P o 

~ 

where p is far-field acoustic mean square pressure extrapolated 

to a radial distance of one meter from the source, the reference pressure 

-5 2 

being 0.0002 microbar (2 x 10 N/M ) for underwater noise measurements. 

If the force exerted by the body on the fluid can be determined then 
the SPL can be calculated. 



CHAPTER III 



HYDRODYNAMICS AND CONTROL SURFACE OSCILLATION 

Principles of incompressible, non- cavitating hydrodynamics 
can be utilized to determine the external forces acting upon the 
fluid which will yield the sound pressure field. Sound generated by 
the oscillating control surface only is considered. This sound pressure 
will be additive to any presently inherent in vehicle operation exclusive 
of the control surface system. 

In Equation (2.2), F^(t) must be determined to find the SPL. This 

total force exerted by the body on the fluid is due primarily to dynamic 

lift directed normal to the approach flow. A secondary force is due to 

viscous drag which is oriented parallel to the approach flow and generally 

equal to about one twentieth of the lift in magnitude. 

For an oscillating hydrofoil the oscillating component of lift is a 

function of approach velocity, hydrofoil geometry, fluid density, and the 

magnitude and frequency of the oscillation. For a steady hydrofoil the 

dynamic lift is a function only of approach velocity, foil geometry, 

fluid density and angle of incidence with approach flow. A frequently 

(5) 

used dimensionless parameter is the lift coefficient C . 

L 

C T = —~ 7 = Kq ( - 17°< a < 17°) 

L £ S u 2 

for a symmetrical thickness form. 
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where L» 

P 

S 

U 

a 

K 



dynamic lift 
fluid density 

area of foil in the span-chord plane 
approach flow velocity 
radian angle of attack 

proportionality constant dependent on aspect ratio 



The aspect ratio is the ratio of span squared to foil area S. For two 
dimensional flow (infinite aspect ratio), K = 2tt. For zero aspect ratio 
(chord length much greater than span), K = ir /2. It is known that at 
reduced frequencies greater than 0. 5 the oscillating component of 
the flow becomes two-dimensional. 

= reduced frequency 

= radian frequency of oscillation 

- semi-chord length 
= approach flow velocity 



K = 



COC. 

TU 



(x> 

c: 

2 



U 



In the course of a mathematical investigation of wing flutter 
phenomena, Theodorsen^^ has shown that 



osc 



J stdy 



= C(k) = 



h ( 2) (k) 



H ( j } (k) + jH^ Z) (k) 



for a lifting surface in rotary oscillation. 

where ^osc ~ ‘complex value of oscillatory lift component 



J stdy 



= quasi-steady lift corresponding to 



maximum inclination 



H 



(2) _ 



n 



Hankel function of second kind of order n 
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The function C(k) is plotted in Fig. 3. 1. Its behavior is quite 
similar to that of a lag network. The reduced frequencies of the 
vehicle control surface in the range of interest were found to be greater 
than 1 . 0 so that C(k) was assumed equal to [0.5 + j 0. 0] throughout 
the analysis . 

The foregoing conclusions have been predicated on the assumption 
that the foil is oscillating in a sinusoidal fashion with time. 

Referring to Fig. 1 . 2(a), assume that the clutch is perfect in that 
the output member locks on with the appropriate input member at the 
switching instant. The angular velocity of the clutch output member will 
be a square wave with half-cycle duration determined by the pulse width 
modulator input. A Fourier decomposition of this signal yields: 



oo 






F(m) e 



2A(l-e 1 ) 

: % > 



m = - oo 

-jmco r 



jmco 2 



m / 0 



F(m ) =< 






(2 -ij, - 1)A ; m = 0 



where 



Q 



1 



c 



clutch output member angular velocity 
switching rate (rad/sec) 



A 



magnitude of input member speed 



T 



duration of positive half-cycle of PWM output 



m 



harmonic index 



T 



period of PWM output 



I 
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Fig- 3-1 Theodorsen's Function 
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In 



1 

the case of steady state zero position error, t = — — 



The 



magnitude of the fundamental component of clutch output velocity will 

, • f 4A 

be a maximum of 



u> 



1 



Referring again to Fig. 1.2(a) relations between control surface 
output velocity and clutch velocity can be written; 






sZ s ^ J 



K , 



where 



Z - load impedance 
n = reduction ratio 
s = complex frequency 



For switching frequencies above load resonance (25 cps ), inertia 



terms dominate and Z — sj . 



Recombining : 



Q 

c 



= n 



1 + 



s 2 (J f + — z ) ^ 

n 

K , 



Examining the relative magnitudes of inertias; 

a 3x10 ^inlbsec^(3-5 J^. , see Fig. 1.6) 
J = 8 x 10 ^ in lb sec^ 

fj 



since a reduction ratio of at least 15:1 will be required, the effects 
of reflected load inertia can be ignored with the result that 






n, 



n [ 1 + 



K, 



indicating that switching frequency harmonics will be attenuated 
inversely as the square of the order. Coupling this result with the — 
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character of the clutch output spectrum, the foil motion can be 
considered a pure sinusoid at fundamental switching frequency. 
The source SPL can now be calculated 



where 



L, 



2 

p 

^ source 



( 



dL ,2 

dt ] 




2 



o 

L sin a>.t = C(k) K T F0 T sintiKtl 
osc 1 ' ' L L L 1 J 

K t = constant of lift-incidence proportionality 

J -i 

= Tr'pS (MKS units) 

C(k) = Theodorsen’s function 



Since the sound pressure is a result of vortex shedding, a 
measure of spanwise phase correlation is necessary. A good measure 
of phase correlation is obtained by a strip theory approach. The foil 
was arbitrarily subdivided into four spanwise sections. The Theodorsen 
function was found to be nearly invariant across the span so that the 
sound field was assumed 100 percent phase correlated. The tip vortex 
baffle increases the reliability of this assumption. 

Converting the mean square pressure to decibel ratio 

SPL = 20 log in (Q t ) + P. dB 

source & 1 0 v L 1 

r q i 



= 20 log 



10 



_ n( 1 -to 



1 K, 



) J 



+ P, 



dB 



P contains K *, C(k), reference pressure (0.002 pb), (0.707) , 
1 L 

16 tt 2 and C 2 . 

o 
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Q 

Q 

If we further assume that = 100 RPM (specified maximum 

n ' 



control surface speed), then 



2 J f 



SPL = P, -20 log. n «,( 1-w. ~~ ) 

source 2 fe 10 P 1 K , ' 



This relationship is plotted in Fig. 3.2 for various ratios of 



switching frequency to 



/ 



K, 
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CHAPTER IV 



CLUTCH POWER DISSIPATION 

As mentioned in Section 1.4 power is dissipated in the interval 
between the switching instant and the point at which clutch input and 
output members become locked. The instantaneous dissipation is the 
product of pullout torque and slip speed. The time average power 
dissipated will be approximately twice the product of clutch output 
inertia, input member speed squared and switching rate. 

In order to minimize power dissipation, the absolute value of slip 
velocity must be minimized. Assuming and switching frequency 

fixed (see Fig. l.Z(s)), the effect on slip velocity of varying was 

observed on an analog simulation. Since will be three to five 
times greater than was assumed infinite for purposes of 

simulation. The system simulated and computing diagram are indicated 
in Figure 4. 1(a). 

The simulation was carried out on the Philbrick Analog Computer 
at the M.I.T. Engineering Projects Laboratory. Results are indicated 
in Figure 4. 1(b). The input member velocity is the square wave. 

The output member velocity is the initially ramp-like square wave 
which tracks the input. The third trace is the amplified difference 
between the two. A small biasing error is evident in seme of the 
computations. This is attributed to amplifier drift since it was necessary 
to saturate one amplifier in the loop. 

From the simulation, it may be concluded that the spring 
should be made as stiff as possible without exceeding clutch pullout 
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COMPUTING DIAGRAM 



Fig. 4.1 (a) Analog Simulation of Clutch Dissipation Phenomena 
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Fig. 4. 1(b) Response of Second Order Clutch Output System 
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torque in spring windup. At first glance, no differences may 
be apparent in Fig. 4 . 1(b). A closer examination reveals that response 

CO CO 

time for ——— =1.21 is about half that for — ■ -• = 0.478. 

CO co 

s w s w 

Input member speed and switching rate must be accounted for 

in the final selection of K^. 
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